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I. INTRODUCTION 



A material point immersed in an environment undergoes a jittering 
movement dubbed Brownian motion. For thermal environments the Brownian 
motion of the tagged particle having position-dependent potential energy 
V(x) may be described by the Fokker — Planck equation in the absence of inertial 
force (the so-called Smoluchowski equation) [1] 



df(x, t) _1 d 
dt ydx 



dV(x) 



dx 



fix, t) 



ksTd^fix, t) 



for the time-dependent probability distribution function fix,t). The diffusion 
constant, given by the Einstein's fluctuation-dissipation relation [2] 

D = ^, (2) 

is expressed in terms of the friction constant y, with dimensions of mass per time, 
and the thermal energy /cgT, where /cg is dubbed the Boltzmann constant and T 
denotes the heat bath's temperature. 

On account of the environmental fluctuations, the escape rate of a Brownian 
particle over a barrier separating two metastable states located at x^ and x^, — in a 
double-well potential, for instance — is said to be a thermal activation 
phenomenon [3 — 9]. On the basis of the Smoluchowski equation (1), Kramers [9] 
found the following steady-state escape rate 

r(oo)=— ksT, (3) 
27zy 

where m is the mass of the Brownian particle, 0)^ the oscillation frequency within 
the bottom well at x^ and (Oi, the oscillation frequency over the barrier located at 

Xi,. 

It is worth noticing that no diffusive motion, D = 0, and no escape rate, 
r(oo) = 0, are predicted at zero temperature. Yet, as far as quantum effects are 
concerned, metastability phenomenon must be characterized as a quantum 
tunneUng process [8,10 — 11]. Accordingly, the diffusion coefficient (2) as well as 
the escape rate (3) cannot vanish at zero-temperature realm. 

Because of the wide practical and theoretical interest in chemical physics, 
for instance, recent studies have been come out showing how quantum effects 
upon the Smoluchowski equation (1) and the Kramers rate (3) can arise. For 
example, in the wake of the Caldeira-Leggett approach [12], Ankerhold etal. [13 — 
16] have derived a quantum Smoluchowski equation and explored its physical 
meaning applying it to chemical reactions, mesoscopic physics, and charge transfer 
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in molecules. Meanwhile, further works [11,17,18] have revealed that such a 
quantum Smoluchowski equation is plagued by some drawbacks, for it may violate 
the Second Law of Thermodynamics, for instance. Another Hamiltonian account 
aiming at to quantize the Smoluchowski equation (1) has been developed 
elsewhere [19,20], whereby a quantum Smoluchowski equation has been derived 
and the problem of quantum tunneling has been studied at zero temperature, in 
contrast to Ankerhold et aJ.'s survey. This approach does not depend on the path 
integral formalism, but it is based on the canonical quantization. Other quantum 
Smoluchowski equations have been also derived in Refs. [21 — 25]. 

Alternatively, we have put forward a quantization method (the so-called 
dynamical quantization [26 — 32]) without alluding to the existence of an 
underlying Hamiltonian model describing the interaction between the tagged 
particle and the environment. So the main purpose of the present paper is to 
explore the physical significance of our method of dynamically quantizing the 
Smoluchowski equation and cope with the problem of tunneling of a Brownian 
particle immersed in a quantum environment without reckoning with certain 
assumptions made in our previous work [32] and criticized by Ankerhold et al. 
[37]. 

Our article is then laid out as follows. In Sect. II we achieve the non- 
Markovian generalization of the classical Smoluchowski equation in terms of the 
concept of diffusion energy, whereas in Sec. Ill we derive a quantum non- 
Markovian Smoluchowski equation for a harmonic oscillator and find out a steady 
solution on the basis of which we investigate quantum tunneling in the long-time 
regime in Sec. IV. At last, summary and perspectives are presented in Sec. V. In 
addition, some technical details are presented in three appendices. 
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II. DERIVING THE NON-MARKOVIAN SMOLUCHOWSKI 
EQUATION 

A. The Einstein — Langevin stochastic approach 

As is well-known, in the classical domain a given isolated system with mass 
m is described by the position q = q(t) and the linear momentum p = p(t) = 
mdq/dt, where the parameter t denotes the evolution time and the symbol d/dt a 
differential operator actuating upon q. Its deterministic dynamics fulfill 
conservative Newton's equations in phase space [33] 

dp dV 

dt^~d^' ^"^^ 

dq p 
dt^m' 

where the inertial force (4) is the position derivative of the potential energy 
V = V(^q) inherent in the particle's motion. 

In order to investigate the erratic motion of a tagged material point 
immersed in an arbitrary environment (a paradigmatic example of open system), 
Langevin [34] extended the Newton's deterministic approach by putting forward a 
random generalization of the Newtonian dynamics (4) and (5) through a set of 
stochastic differential equations (the so-called Langevin equations ) 

dP dV(X) P 

-7- + 7mt), (6) 



dt dX m 

dX _P 
dt m 



(7) 



where the environmental force. 



F,^^{p,^') = -Y^ + Yb^'{t), (8) 

is made up by a dissipative force —yPlm — exhibiting no memory effect and 
accounting for stopping the particle's motion — and by a fluctuating force, yhWit), 
which is responsible for activating the particle's movement. The parameter 
h controls such an environmental influence (fluctuations), while the 
phenomenological parameter y, which is expressed in dimensions of mass per 
time, is termed the damping constant. The term hW{t) in (6) is in dimensions of 
velocity, i.e., [length x time~^]. Thus, we can readily check that the parameter b 
may be expressed in dimensions of [length x time~^^^] and the function "/^(t) in 
dimensions of [time"^/^]. The existence of the environmental force (8) suggests a 



4 



relationship between the dissipating and fluctuating processes, i.e., between b and 
7, since < 7 < oo and < 6 < 00. 



From the mathematical viewpoint the quantities X — X(t), P — Pit), and 
V = V^(t) in the Langevin equations (6) and (7) are viewed as random variables in 
the sense that there exists a probability distribution function, Txpiy(.x,p,il},t), 
associated with the stochastic system {X, P, V], expressed in terms of the possible 
values X = [Xi], p = {pj, and ifj = {ifji}, with i > 1, distributed about the sharp 
values q, p and (p of X, P, and V, respectively. 

In the Einstein — Langevin randomization process the average values of X, P 
and V are expressed as 

+ CO 

(X) — jjj xTxpw (x, p, xp, t)dxdpdip, (9) 

— CO 

+00 

<P) = IW pTxpw(.x,p,xp,t)dxdpdxp, (10) 



+00 

(W) ^ jjj ipTxpw(x,p,ip,t)dxdpdip. (11) 

—00 

In the deterministic case, i.e., as b = and J^xpvCxiP.tp.t) = S{x — q)d(p — 
p)S(jp — q)), the deterministic variables x = q and p = p = mdx/dt are recovered 
as the respective average values (X) — q and (P) = p. 

1. The generalized Smoluchowski equation 

We start from the Langevin equations (6) and (7) in the Smoluchowski limit 
(also known as the large friction case) 



\FenviP.'y)\ » 



dP 



dt 



(12) 



such that inertial effects in (6) can be negligible, i.e., dP/dt = 0. So the focus on the 
motion of the Brownian particle is gained through a single stochastic differential 
equation given by 

dX dViX) 

in which the dissipative force ydX/dt is offset by the conservative force, 
-dV{X)/dX, and the stochastic force, ybW^t). 
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Noticing that 



t+E 



lim f {W(t'))dt' = 0, (14) 

£^0 J 
t 

the Langevin equation (13), divided by y> gives rise to the following Fokker — 
Planck equation on configuration space in the Gaussian approximation (see 
Appendix A) 



dfix, t) _ d 
dt dx 



1 dVix) 



d f(x t) 

f{x,t)+'D{t)-j^. (15) 



where the mean value of 'F(t) is given by 



t+E 



<nt)) = lim-[ <nt'))dt', (16) 

t 

and the time-dependent diffusion coefficient 

2)(t) = -y- (17) 
is expressed in terms of the following function 

2 2 

S{t) = ^/(t) = ^lim- ff {Vit'mt"))dt'dt", (18) 

Z L £^0 £ jj 



which has dimensions of energy, i.e., \mass x length^ x time~^]. Hence we call 
£(t) the diffusion energy responsible for the Brownian motion of the particle 
dipped in a generic environment. 

The time-dependent function /(t) in (18), which encapsulates all the 
fluctuation effects on the Brownian particle through the auto-correlation function 
{b^'(t')bW(t")), is a dimensionless term because b and ^(t) in (6) has been 
assumed to have dimensions oi{length x time~^^^] and [time~^^^], respectively. 

An outstanding feature underlying our diffusion energy concept in (17) is 
that £(t) = 7'D(t) conveys in tandem the fluctuation phenomenon through the 
diffusion coefficient D(t) as well as the dissipation process encapsulated in the 
friction constant y. We can then state that our definition of diffusion energy (17), 
which fulfils the validity condition 

< £(t) < oo, (19) 

stands out a general fluctuation — dissipation relation underlying all open systems 
described by the Langevin equation (13) and its corresponding Fokker — Planck 
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equation (15). Both cases £(t) = and £(t) = oo are not concerned, for they may 
violate the fluctuation — dissipation relation. The former case may lead to 
dissipation without fluctuation, while the latter one may give rise to fluctuation 
without dissipation. 

In the long-time regime we assume that the diffusion energy (18) renders 
steady, i.e.. 



with 



Um £(t) = £(cx)) 

t^oo 



lim /(t) = /(oo) = 1. 

t-»oo 



(20) 



(21) 



The physical significance of condition (21) has to do with the fact that the 
environmental fluctuations do possess Markovian correlations, i.e., /(t) displays a 
local (short) range behavior decaying to one in the steady regime. In generic open 
systems non-Markovian effects show up at the range < t < oo. The function /(t) 
exhibiting the asymptotic behavior (21) can be built up as (see Appendix B) 



/(t) = 1 - e' 



(22) 



such that at t = the diffusion energy associated with the initial probability 
distribution function f{x, t = 0) = 5(x) vanishes: £(0) = 0. In (22), denotes 
the correlation time making the fluctuations non-Markovian. 



Accordingly, from (18) we obtain the fluctuation-dissipation relation in the 



form 



b = 



2£(oo) 



N 

and the Fokker — Planck equation (15) reads 

dVix, t) 



a/(x, t) _ 1 5 
dt Y dx 



dx 



■f{x,t) 



dx^ 



with the effective potential 

V{x, t) = V{x) - Xyl2Ye(oo){W(t)). 



(23) 



(24) 



(25) 



For the case in which {^(t)) = and /(t) = 1, the Fokker — Planck equation 
(24) reduces to the Markovian Smoluchowski equation [1] 



a/(x, t) _ 1 a 

dt y dx 



dVix) 



dx 



f{x,t) 



+ 



£(cx))aV(x,t) 

V dx^ 



(26) 
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Thereby, we dub Eq. (24) the generahzed Smoluchowski equation because it 
reckons with both non-Markovian and averaging effects. 

III. DYNAMICAL QUANTIZATION OF OUR NON-MARKOVIAN 
SMOLUCHOWSKI EQUATION 

A. Quantum non-thermal open systems 

Our classical Smoluchowski equation (24) with (22) at point x for the case 
of a harmonic oscillator, V{x') = /cx'^/2, where the constant k has dimensions of 
[mass X time~^], reads 



dfCx, t) _k d 
dt Y dx 

in terms of the position 



[xfix,t)]+—^[l-e , (27) 



x = x'-^V2rf(oo)<nt)>. (28) 

Our purpose is to quantize the non-Markovian Smoluchowski equation (27). To 
this end, we start with it at points Xi and X2, i.e., 

aj(xi,t) k kx^dxix^.t) £{t)d^xi^i,t) 

= -Xixi, t) + — — + — ^ (29) 



dt Y Y dxi Y dxl 



and 



dx{x2,t) k kx2dx{x2,t) £{t) d'^xi^i.t) 

^7 = -Xix2, t) + + . (30) 

dt Y Y 9x2 Y 9x2 



Solutions xixi,t) and xix2,t) are deemed to be associated with the diffusion 
energy £(t) = 2f(cx))(l - e"^/^0- By multiplying (29) by x(x2,t) and (30) by 
x(.Xi, t) and adding the resulting equations, we obtain 

where f = f (xi,X2, t) = xi^i, t)x{x2, t) = V/(^i< 0- 

By performing the following change of variables into configuration space, 
(xi,X2) ^ (x.ryh), 

77h 

Xi = X - y, (32) 

jyh 

X2 = X + y, (33) 
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and making use of the relations 



d Id Id 



d Id Id 



dxi 2dx hdt] 



dx2 2dx^ h drj' 



the classical equation (31) changes into the quantum equation of motion 



dp 2k k dp k dp £^(1) (d^p 4 d^p 




(34) 



dt Y Y dx Y dri y \dx^ drp- 



where we have replaced the solution ^ = ^{xi,X2,t) with p = p{x,r],t) and 
employed the notation 



to display the quantum nature of the diffusion energy which turns out to be now 
expressed in terms of the Planck constant h. Notice that the constants k, y, and the 
dimensionless function /(t) hold classical during the classical-quantum transition 
(32) and (33). 

The geometric meaning of our quantization conditions (32) and (33) is 
related to the existence of a minimal distance between and X2 due to the 
quantum nature of space, X2 — Xi — r/h, such that in the classical limit h ^ 0, 
physically interpreted as 



the result X2—Xi=x can be recovered. Accordingly, the classical limit of the 
quantum equation of motion (34) leads rightly to our Smoluchowski equation (27) 
because 




(35) 



|77h| « \X2 -Xil, 



(36) 




(37) 



provided that 



lim£h(t) = £(t). 



(38) 



By making use of the Fourier transform 



00 




(39) 



—00 
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which changes the variables from quantum configuration space (x,r]h) onto 
quantum phase space (a:, p; h), equation of motion (33) turns out to be written 
down as 



dW _k dW k dW £t,(,t)d^W 
dt dx dp y dx^ 



4£h(t) 2 



y h^y 



W (40) 



in terms of the function W = I/K(x, p, t). Because the exponential factor e''"' in (39) 
is to be a dimensionless term, it follows that the variable p is to have dimensions of 
linear momentum. 

On account of the following property (see Appendix C) 

lim Wix, p, t) = T{x, p, t) = fix, i)5{p), (41) 

the classical limit of (40) leads to the Smoluchowski equation (27)as well. 

1. The initial condition 

Our quantum phase-space Smoluchowski equation (40) may be solved 
starting from the non-thermal initial condition 

Wix.p, t = 0) = — e V h (42) 
Trn 

leading to the distribution F{x,p) = S{x)S{p) in the classical limit, h ^ 0. The 
constant a has dimensions of time per mass, that is, we may take a = l/mo), 
where o) = -Jm/k is the oscillation frequency of the harmonic oscillator. The 
probability distribution function (42) generates the fluctuations 

T hmo) 

(P^) = (43) 

which are compatible with the Heisenberg indeterminacy principle 

WW^) = I- (45) 

The initial mean quantum-mechanical energy of a Brownian harmonic oscillator is 
given by 

albeit its initial diffusion energy (35) is null, i.e., £^(,0) = 0. 
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2. The time-dependent solution 

A time-dependent solution to (40) reads 



1 r i ^ghCQp^ i 

W(x,p,t) ^—e i4Bh(t) J, (47) 

Trn 



where the function 



_2fci r h /I M 1 t^e tc \ 



4ma) ''^ ^V2fc (y - 2fctc)/J \2k (y - 



is expressed in terms of the evolution time t, the relaxation time tj. = y/k, and the 
correlation time responsible for non-Markovian effects. 



3. The long-time regime 

In the long-time regime, t ^ oo, physically interpreted as t » t^, t^., the 
probability distribution function (47) renders steady, that is, 

1 r fex^ 2£n(cx.)p2l 

l^stUp) = -re ^2^'^^°°) fch' J (49) 
Trh 



leading to the fluctuations 



/ch2 

a^> = ^, (51) 

compatible with the Heisenberg relation (45), and the steady mean quantum- 
mechanical energy 

which is to be defined for the quantum diffusion energy satisfying the validity 
range < S^ipo) < 0. In the especial cases S^ipo) = and £^(,oo) = oo expression 
(52) blows up, thereby violating the quantum dissipation-fluctuation relation. 
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B. Quantum thermal open sj^ems 

Imagine that the environment is made up of a set of harmonic oscillators 
having the same oscillation frequency o), so that the mean energy of this quantum 
heat bath after attaining the thermodynamic equilibrium at temperature T is given 
by [35] 



E = — — r = — coth 

,ekBT -1/ « 

where the h-dependent energy, a)h/2, corresponds to the mean quantum- 
mechanical energy (46) of the heat bath prepared at the initial state (42) and 
k^T the classical thermal energy of the heat bath, where kg is the Boltzmann 
constant. 

We now assume that the Brownian particle's quantum diffusion energy 
£ji(oo) can be identified with the heat bath's quantum thermal energy E , i.e., 



, ^ coh / coh \ 
£,(oo) = _coth(— ). (54) 



(i)h _ ( (xih 

Using (35) in the steady regime the quantum diffusion coefficient can be derived as 

B,(cc)=-coth(— ). (55) 

The steady probability distribution function (49) becomes the steady 
quantum-mechanical equilibrium function 



WAx,p) = ^e'-''^w) (56) 



leading to the fluctuations (50) and (51) in the form 



, mo)h 

{P ) = (57) 



h _ / a)h 

The mean quantum-mechanical energy (52) is then given by 

(oh ojh _ ( (lih 



h ( (jdh \ 



a)h (xih ( (lih \ 

(E(cx))) = ^+ coth —— . (59) 
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1. Zero temperature 

The steady probability distribution function (56) at zero temperature reads 



W,Y=°\x, p) = — e'y^^i (60) 



and the fluctuations (57) and (58) become 



(^-2^(7=0) ^ (62) 

The mean quantum-mechanical energy (59) is then given by 

<E(cx)))(^=o) = ^ (63) 

which equals the initial energy (46) and the diffusion energy (54) in the steady 
regime 

£i^=°>(oo)=f, (64) 
whereas the quantum diffusion coefficient (55) turns out to be 

j'r°'(»)=f- (65) 

On account of the energy conservation (54), results (60 — 65) provide 
physical significance to our quantum phase-space Smoluchowski equation (40) at 
zero temperature. 

2. The classical limit 

The classical limit h ^ should be physically interpreted as the classical thermal 
energy is deemed to be too large in comparison to the quantal energy, i.e., 

a)h 

k^T » (66) 

Thus, our quantum Smoluchowski equation (40) reproduces the classical 
Smoluchowski equation (27) at the high temperature range, whereas the classical 
limit of (59) leads rightly to the principle of equipartition of energy 

lim<E(cx))) = = {Eioo)), (67) 
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while the quantum diffusion energy (54) turns out to be the thermal energy 



£(cx)) = keT. (68) 
We also obtain the Einstein's celebrated relation 

2)(cx)) = (69) 

as quantum effects can be negligible in the diffusion constant (55). In other words, 
the Einstein diffusion constant (69) maintains valid at high temperature 
T » a)h/2/cB. 

Furthermore, we can derive the Maxwell — Boltzmann distribution function 

fix) = p~W (70) 

^InkBT 

as the classical limit of the (marginal) probability distribution function (56). 

In view of the energy conservation (54), results (66 — 70) yield physical 
meaning to our quantum phase-space Smoluchowski equation (40) at high 
temperature by reproducing the well-known tenets of the classical theory of non- 
inertial Brownian motion at thermodynamic equilibrium. 



IV. AN APPLICATION: THE QUANTUM KRAMERS EQUATION 

A. Quantum non-thermal open systems 

In order to calculate the quantum Kramers escape rate of a Brownian 
particle over a potential barrier in the strong friction regime, we consider as a 
starting point the steady solution (49), i.e., 

1 r i^(x) 2£ft(cx,)p2 l 

W,^ix,v)=-re ifihC"')^ feh^ J. (71) 
nn 

We perform an expansion of the potential function V{x) — in a Taylor series 

around a certain point Xq, i.e., V{x)~V{xq) + {k/2'){x — Xq)^, and construct the 
IVlarkovian non-equilibrium function 

[y(xo)+|(x-xo)^] 2£a(°°) 2 
VK>eq U P) = «<P(^- P)e ^h^"") , (72) 

a being a constant and cp = (p{x,p) a function to be determined. On inserting 
solution (72) into our quantum Smoluchowski equation (40), given by 
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-(^-^o)^ P p 

y ox y op 



neq 



k 4£h(oo) 2 
7 7h2 ^ 



M^neq = 0, (73) 



with M^eq — l^eq p)> we obtain the equation of motion 



'neq 



k d(p k d(p £f^(^oo)d^(p 2k dcp 

y ox y op y ax^ y ax 



(74) 



which changes into 



dcp 



in terms of the variable 



^ = P(x-Xo)-p. 



(75) 



(76) 



The constant p in (75) and (76) has dimensions of [mass x time~^ ] such that the 
term P(x — Xq) in (76) is in dimensions of momentum, Hke p and f . A solution to 
the differential equation (75) reads 



(77) 



.27r)S2£^(oo) 

—00 

fulfilling the condition ^(f ^ 00) = 1. This result (77) is only possible if the 
potential curvature is negative, i.e., /c < 0. Hence, we hereafter use Xq = x^,, and 
k = —moil, where the quantity denotes the particle's oscillation frequency over 
the potential barrier at x^,. So the barrier top is located at point x^, while the two 
bottom wells are at x^ and x^, such that Vipc^ ) = = V{Xa), with x^ < X/,. 



Using (77), the stationary function (72) becomes 



VKieq(^<P) = « 



2£h(oo)„2 



from which we can find out the probability current as 

CO 

acolh^m 



/ 



e Wsr^d^ (78) 



h 



m J 



P^neqix = Xb.p) dp = - 



_VOcbl 



4£h(oo) ^m^h^a)t-4p^el{oo)' 



(79) 



where we used the result 



00 ?=p 

f 2 r f2 1 / TT 
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u and V being non-negative constants. 

The number of particles through the well located at is calculated as 



ria = U Weq(x,p)dpdx = anh, 



(80) 



whereby we used the equilibrium distribution function 

coming from the non-equilibrum function (72) for (p(x,p) — 1 calculated at the 
well Xa, with ka — mcol. 

Making use of (79) and (80) the quantum Kramers escape rate of a 
Brownian particle immersed in a non-thermal environment is then written down 
as 

r(oo) = ^ = — ^ ^ (82) 

in terms of the steady quantum diffusion energy, £h(°°) — 7'^h/2 — 7l'h(°°)' that 
fulfills the condition 

< £h(oo) < (83) 

It is worth stressing that our main upshot (82) does display dimensions of inverse 
of time, as it should be expected, as well as being compatible with the dissipation- 
fluctuation relation S^ipo). 

Our quantum escape rate (82) decays exponentially with respect to the 
ratio— V(x/,)/£ji(oo). In addition, the so-called Arrhenius factor, given by 

cofh^m 1 
a - — ^^^^^=^^^^= (84) 

does depend upon the arbitrary parameter p introduced by the change of variables 
(76). Because it is expressed in dimensions of mass per time, we may consider the 
following cases: P = Y and P = Nmo)},, where TV is a dimensionless number such 
that < yv < 00. 
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1. TheP = Nmoif, case 

Considering (i — Nmcoi, our quantum rate (82) turns out to be given by 



r(oo) = 



V(.Xb) 



with 

Taking TV = in (85) leads to 

r(oo) = — -V^e 

47r£h(oo) 

2. 7»e/? = Y case 

For P =Y, our rate (82) reads 

r(cx)) = 



4n:£h(oo) ^Tn2h2a)^-472f2(cx)) 
provided that 

< £h(cx)) < ^— . 
B. Quantum thermal open systems 

For thermal systems the quantum diffusion energy S^ipo) given by 
evaluated over the potential barrier, i.e., 



£h(oo) = ^coth(^j. 



Accordingly, our quantum Kramers rate (82) reads 



a,jhcoth(ig^) 



2- (^) Jm2a,|-;?2,oth2(i^)" 
valid for the low-temperature regime 



0<coth -A:- < 



At zero temperature, Equation (90) becomes £hi°°) = o)h/2 and (91) reduces to 

_2V(xb) 

,„ , mcol e 
r(^=o)(cx)) = ^ V =■ (92) 

27r ^m^col-p^ 



1. Thep = Nmcdi, case 

For = Nmco\), we find from (91) the following quantum escape rate 

__m^i^)__ 

oifthcoth ( ^r^l 

r(^) = ' , , . _ (93) 



with 



< yv < ^ ■ (94) 

The N = case provides 

r(oo) = e (95) 

2;rC0th(^) 

At zero temperature, (93) reduces to 

_2V0cb) 

rO'=o)Coo) = ^ _ (96) 
with < yv < 1. As yv = we have 



r(^=o)(oo) = -^e ^fft . (97) 
2ti 
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2. Thep = Y case 



For P — y, Eq. (91) yields 



2- (^) y^l_y2,,,^2(^^'^ 



r(^) = r\.^^ ^ (98) 



and 

2V{xh) 



r(^=o)(cx)) = V (99) 

A remarkable feature underlying our results is the absence of any damping 
constant y ii^ the quantum Kramers rates (93), (95 — 97). This fact thereby 
suggests the existence of a quantum transport phenomenon coming about without 
dissipation (superfluidity) at thermal equiUbrium without violating the quantum 
dissipation-fluctuation relation (54). This counter-intuitive quantum tunneling 
process in the strong friction regime shows up as a consequence of the arbitrary 
character of the parameter p in the linear transformation (76) performed on our 
quantum Smoluchowski equation (74). That result is the same as the one obtained 
in Ref. [32] and has been already reported in studies about the influence of 
dissipation on the Landau — Zener transition at zero temperature [36]. 

It is worth pointing out that our non-dissipative quantum Kramers rate (93) 
can be defined at 7 = in accordance with our previous work [32] and the Indian 
group's approach [19,20], but in contrast to Ankerhold's approach [13 — 16]. 

As far as the P = y case is concerned, our quantum tunneUng escape rate 
(98) and (99) can depend upon the friction constant, thereby predicting no 
superfluidity phenomenon in the large friction realm. This feature seems to be in 
accordance with the works in Refs. [19,20, 37]. 

C. The breakdown of the dissipation-fluctuation relation 

Although our quantum Kramers escape rate (82) cannot be defined for £^(,oo) = 
without violating the dissipation-fluctuation relation, we may conjecture about the 
mathematical behavior of (82) as £^(,co) corresponding to the extreme 
physical situation 



8m£ji(oo) 



» 



(100) 



Under this condition the mean quantum-mechanical energy (52) approximately 
equals the mean kinetic energy, i.e., 
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{E(oo)) - 



8m£ji(oo)' 



(101) 



From (82) in the limit £h(°°) ^ 0, we obtain then the result 



r(oo) = cSiXh). 



(102) 



The constant c having dimensions such that r(oo) holds the dimensions of inverse 
of time. The non-exponential behavior of the escape rate (102) may suggest a 
breakdown of the dissipation-fluctuation relation provided that the condition 
(100) is vahd, i.e, £ii(oo) « o),,h/2. 

Experimental evidence for tunneling processes undergoing non-exponential 
decay have been reported in the literature in the context of isolated quantum 
systems [42]. Yet, our theoretical prediction (102) has to do with quantum open 
systems in the strong friction regime. So, we hope our outcome (102) could 
prompt experimental researches to investigate quantum tunneling processes in 
the presence of dissipative open systems in the large friction regime. 



In Sec. II we have derived the non-Markovian Smoluchowski equation (24) 
expressed in terms of the fluctuation-dissipation relation, given by the diffusion 
energy (18), and the averaging effects giving rise to the effective potential (25). 

In Sec. Ill the dynamical quantization of our generalized Smoluchowski 
equation (27) for the case of a Brownian harmonic oscillator has been achieved. In 
so doing, our quantum phase-space Smoluchowski equation has been obtained 
[Eq.(40)] in terms of the quantum diffusion energy (35). In the long-time regime 
the quantum non-thermal stationary solution (49) has been derived. 

As far as a quantum thermal bath is concerned, the quantum diffusion 
coefficient (55) associated with the quantum diffusion energy (54) has been found 
out. We have also obtained the steady thermal quantum probability distribution 
function (56). All our results are valid at zero temperature, 7 = 0. On the other 
hand, the Einstein's renowned diffusion constant (69), the energy equipartition 
(67), and the Maxwell — Boltzmann probability distribution function (70) have 
been rightly derived in the classical limit. 

In Sec. IV the non-equilibrium quantum Kramers rate (82) for steady non- 
thermal open systems has been found on the basis of our quantum Smoluchowski 
equation (40). Our upshot (82) is expressed in terms of the dissipation-fluctuation 
relation, which is given by the steady quantum diffusion energy £^(,co) = 7!Dji(oo), 
as well as in terms of the arbitrary parameter p introduced by the linear 
transformation (76). We have also predicted a non-exponential behavior of our 
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quantum escape rate (82) as £^(00) 0, i.e., on the edge of the breakdown of the 
dissipation-fluctuation relation. 

For thermal quantum systems our main finding is that the tunneling rate 
(88) may depend upon no friction constant [Eqs. (93,95 — 97)]. This fact hereby 
suggests the existence of a damping-independent tunneling process in the strong 
friction regime. Yet, supposing )S = 7 in (91) leads to (98) and (99) in which the 
quantum tunneling phenomenon depends on the damping constant. 

Lastly, it would be quite catchy to compare our quantum phase-space 
Smoluchowski equation (40) and its predictions with other accounts existing in 
the literature [13 — 25] and bring up new features underlying our non-Hamiltonian 
approach to open systems by applying it to the Brownian motion of particles 
dipped in fermionic and bosonic thermal reservoirs on the ground of our concept 
of quantum diffusion energy £^(,co). This task will be achieved in a forthcoming 
paper [38]. 
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APPENDIX A: GENERALIZING THE SMOLUCHOWSKI 
EQUATION 

The Langevin equation (13), given by 



dX IdV 

-J-= -j7;+b^it), 

dt ydX 



(Al) 



gives rise to the Kolmogorov stochastic equation on configuration space [39 — 41] 

df(x, t) 



dt 



= Kf{x,t), 



(A2) 



where the Kolmogorovian operator K acts upon the function f(x, t) according to 

on , , 



Kf(x, t) = ^ ^ [Ak (x, t)/(x, t)] , 

fc=i 



(A3) 



the coefficients Aj^ (x, t) being given by 

mx)^) 



A^ix.t) = lim 

E->0 



(A4) 



where the average values, ((AX)''), in the coefficients (A4) are to be calculated 
about the sharp values q, i .e.. 



Txw ix, \p,t) = S{x- q)T,p {\p, t). 



(A5) 



According to the Pawula theorem [40], if the coefficients Aj^i^x, t) in (A4) are 
finite for every k and if Aj^(x, t) = for some even k, then Aj^^x, t) = for all k > 
3, thereby assuring the positivity of the probability density function f(x, t). So, if 
(a:2 — a:i)^~0 such that (^^2—^:1)^ = 0, it follows then that A]f{x,t) = 0, k> 
3, where X2 = x{t + e) and Xi = x(t). Accordingly, the Kolmogorov equation (A2) 
can be reduced to the generalized Smoluchowski equation (the so-called Fokker — 
Planck equation in configuration space) 



df(x,t) d 1 

-^^g^ = -g^ [Ai(x, t)f{x, 0] + 2 ^ [^2(^, t)f{x, t)l 



(A6) 



where the drift coefficient is given by 



Al {x, t) — lim 
£^0 



IdV 



(A7) 



and the diffusion coefficient reads 
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A2 {x, t) = lim 

£^0 



t+E 

= ^^lim jjj {Vit')^it"))dt'dt", 
t 



(A8) 



with 



t+s 

lim <nt'))dt'=<no) 



(A9) 



and 



t+s 

lim f {1'{t'))dt' = 0. 

t 

To evaluate the coefficients (A7) and (A8) we have used (Al) in the form 

t+E 

edV r 

AZ = Z(t + £) - Z(t) = — -;T7p+b\ nt')dt'. 

Y dX J 
t 



(AlO) 



(All) 
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APPENDIX B: THE TIME-DEPENDENT DIFFUSION ENERGY 



The time-dependent diffusion energy in (18) 

S(t) = £(oo)/(t) (Bl) 
presents the correlational function 



t+£ 



/(t) = lim- ff {Vit')V(it"))dt'dt" . (B2) 

e^O £ JJ 



t 



On the condition that the autocorrelation function, {W(t:')W(t")), can be 

given by 

/ it'+t")\ 



{V(it')V(it")) = yi-e 2t, js(^t'-t"), (83) 

where is the correlation time of V(t) at times t' and t", it follows that (82) 
becomes 

/(t) = 1 - e tc. (84) 

Non-Markovian effects show up in the dimensionless function (84) in view 
of the presence of the correlation time t^, such that in the steady regime, t ^ oo, 
physically interpreted as t » t^, it reduces to /(oo) = 1, thus making the stochastic 
process Markovian. 
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APPENDIX C: PROPERTIES OF THE QUANTUM FUNCTION (39) 

The quantum function W(x,p,t) defined in Eq. (39) does possess the 
following properties: 

(i) just as p{x, T], t) the function W{x, p, t) is deemed to be real; 

(ii) the classical limit provides 



lim W{x, p, t) = fix, t)5{p) = T{x, p, t), 



(iii) the normalization condition reads 

00 00 00 

\\wix.v.t)dxdv= \\pix.n.t)mdnd.= \pix.t)dx = l: 

— 00 — 00 — 00 

(iv) in the classical limit the marginal distribution 



jwix.p.mp=mx.t) 



leads to 



lim W(x,t) = f{x,t); 

(v) in quantum phase space the Fourier transformation (39) reinstates the 
stochastic character of the momentum p, i.e., 



\w(x,v.t)dx=wxv.ty 

— 00 

(vi) In the classical limit, h ^ 0, we find the result 



Um W'{p, t) = lim | W(,x,p, t)dx = Sip), 



meaning that in the classical domain the momentum turns out to be a deterministic 
variable in the large friction regime. 
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On the basis of the properties above, we could mathematically interpret the 
function W(^x,p,t) in (39) as a quasi-probability density, for it may assume 
negative values in some regions of quantum phase space (x,p;h). Even so, the 
average value of any quantum physical quantity <A(X,P,t) can be deemed to be 
calculated as 



00 



{Jl{X,P,t))^ cA{x,p,i)W{x,p,t)dxdp. 




— 00 
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